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Argelander-Institut fiir Astronomie, Universitat Bonn, Auf dem Hiigel 71, 53121 Bonn Germany 

ABSTRACT 

Context. The Cepheid mass discrepancy, the difference between masses predicted from stellar evolution and stellar pulsation calcula- 
tions, is a challenge for the understanding of stellar astrophysics. Recent models of the eclipsing binary Cepheid OGLE-LMC-CEP- 
0227 have suggested that the discrepancy may be resolved. 

Aims. We explore for what physical parameters do stellar evolution models agree with the measured properties of OGLE-LMC- 
CEP0227 and compare to canonical stellar evolution models assuming no convective core overshooting. 

Methods. We construct state-of-the-art stellar evolution models for varying mass, metallicity, and convective core overshooting and 
compare the stellar evolution predictions with the observed properties. 

Results. The observed mass, effective temperature, and radius of the two stars in the binary system are well fit by numerous combi- 
nations of physical parameters, suggesting a Cepheid mass discrepancy of 10-20% relative to canonical stellar evolution models. 
Conclusions. The properties of the observed binary Cepheid suggest that the Cepheid mass discrepancy is still a challenge and re- 
quires more specific observations, such as the rate of period change, to better constrain and understand the necessary physics for 
stellar evolution models to resolve the discrepancy. 
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1. Introduction 

Classical Cepheids are ideal laboratories for stellar astrophysics 
and evolution. Stellar pulsation provides a measure of the fun- 
damental properties of these stars including luminosity, radius, 
and mass that can be compared to predictions based on stellar 
evolution models. Cepheid pulsation masses, however, d o often 
not agree with predictions from stellar evolution models (IStobid 
11969; Cox 1980; Keller 2008). 

This Cepheid mass discrepancy, defined through the rela- 
tive difference between evolutionary and pulsation masses, is 
currently measured to be about 17% jKeller 2008) and may 
be a function of mass ('Caputo et al.l 120051) and metallicity 
([Keller & Wood 2006). There have been a number of solutions 
suggested to solve this mass discrepancy including convective 
core overshooting in the main sequence progenitors of Cepheids 
and mass loss during the Ce pheid stage of evolution dBono et al.l 
l2006l) . iNeilson et al.l ([201 1) showed that the mass discrepancy 
may be solved by the combination of pulsation-driven mass loss 
during the Cepheid stage of evolution and moderate convective 
core overshooting. All of these suggestions imply the origin of 
the mass discrepancy lies within the stellar evolution models. 

Dynamically measured masses from Cepheid binaries have 
the potential to constrain and verify stellar pulsation or stel- 
lar evolution masses. Dynamical masse s have been determined 
for a ha ndful of Galactic Cepheids ('Evan s & BoltonI Il990t 
lEvans e t al. 1997, 2006, 2008), suggesting that dynamic masses 
are consistent with pulsation masses. The r ecent analysis of 
the echpsing binary OGLE-LMC-CEP0227 bv lPietrzvriski et alJ 
(l20Toh . however, obtains a Cepheid mass of M = 4.14+0.05 Mq, 
the most precise Cepheid measurement to date. The authors also 
derive the stellar pulsation mass for this object to be Mp - 
3.98 ± 0.28 Mq, consistent with the dynamic mass. The authors 
argued that because the pulsation mass agreed with the dynamic 



mass then the Cepheid mass discrepancy can only be solved by 
new physics in the stellar evolution models. However, they did 
not test this assertion with comparisons to stellar evolution mod- 
els. 



ICassisi & SalarisI (1201 ih considered stellar evolution mod- 
els with moderate convective core overshooting and found that 
the models fit the observed effective temperature and radius for 
stellar masses consistent with dynamic masses. As a result, the 
authors conclude that there is no Cepheid mass discrepancy, at 
least for this particular Cepheid. It is debatable whether this is a 
reasonable conclusion. On one hand there is significant evidence 
suggesting the ne ed to include convective overshooting in evolu- 
tion m odels (e.g. ISandberg Lacv et alJl20Tol : iLovekin & Goupill 
2010^. 

On the other hand, the underlying physics of convective core 
overshooting is not understood, hence the physics must be pa- 
rameterized. Furthermore, because overshooting is parameter- 
ized in evolution models, it is not necessarily treated in the 
same way by various evolution codes. Another issue is that 
various observations constraining overshooting seem to r equire 
differing amounts of overshooting (e.g. [K eller & Wood 120061 : 
ILovekin & Go upil'201()':'Sandberg La cy et alj20ia) . Finally, for 
Cepheids in particular, parameterizations o f overshooti ng ma 
act to hide different physics. For instance, [Bono et aO (1200^ 
noted that the Cepheid mass discrepancy might be resolved 
by rotational mixing, mass loss, unknown opacities as well as 
overshooting. These different physical processes are degenerate 
when predicting stellar masses. 

iKelled (l2008h and iNeilson et al] (1201 lb predict a similar 
amount of Cepheid mass discrepancy by assuming different 
physical processes. iKelleil (2008) di d so by assuming only 
convective core overshooting, while INeilson et al.l ( 1201 1 ) pro- 
posed a combination of pulsation-driven mass loss and moder- 
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ate co nvective core overshooting. It should be noted that lKelled 
(I2008h used a convective core overshooting parameterization 
th at is a fact or of two greater than the parameterization used 
bv lNeilson et al. (201 I Q, which is described in the next section. 
In fact, we (iNeilson et al."201 1") misinterpreted the parameteri- 
zation used by Keller (2008), however, this does not affect the 
result that pulsation-driven mass loss can account for 5 - 10% of 
the mass discrepancy, nor the amount of convective core over- 
shooting in that paper. The large range of measured mass dis- 
crepancies and uncertainties determined by ' Keller! (l2008l) are 
still consistent with this result. As such, it is useful to define the 
Cepheid mass discrepancy as the difference between predicted 
masses from stellar pulsation modeling and predicted masses 
from stellar evolution models, assuming no convective core over- 
shooting or other extra mixing process or enhanced mass loss. 
These stellar evolution models can be referred to as canonical 
stellar evolution models. In this way, we can predict a mass dif- 
ference and then use extra physical processes or parameteriza- 
tions to predict a r esolution to the m a ss disc repancy. 

In their work, ICassisi & SalarisI dloTlh made two assump- 
tions about the Cepheid binary: 1) that the Cepheid is evolv- 
ing along the blue loop, and 2) that the Cepheid's metallicity is 
Z - 0.008, with dy/dZ = 1.4. The purpose of this work is to test 
how the predicted stellar evolution masses depend on these two 
assumptions and what this suggests about the Cepheid mass dis- 
crepan cy. As such, we can verify the results of Cassisi & Salarii 
and test the uniqueness of those results. 

2. Stellar evolution models 

We compute model stellar evolution tracks using the 
lYoon&Lange3(l2005h stellar evolution code. The code includes 
prescriptions for mass loss (for hot stars the Kudritzki et alj 
(1989) prescription is used while the deJageretal. (1983) 
prescription is used for cool stars) and convective core over- 
shooting. Overshooting is treated by assuming that for an 
evolutionary timescale convective eddies penetrate some frac- 
tion of a pressure scale height above the convective core, defined 
as dc = UcHp. The coefficient Oc is a free input parameter for 
the evolution code, and is line arly correlate d with the amount 
of Cepheid mass discrepancy (iKeUerll2008h . If one considers 
a Cepheid with a given luminosity then it also has a specific 
helium core mass. The helium core mass is determined by a 
combination of the initial mass of the star and the parameterized 
amount of convective core overshooting. Therefore, if one fits 
the given luminosity with a canonical stellar model, one would 
predict a canonical stellar mass. If the amount of convective 
overshoot is increased, then a smaller initial mass is required to 
fit that given luminosity. Furthermore, one can fit the amount of 
convective core overshooting to fit a predicted stellar pulsation 
mass, meaning that the amount of convective core overshooting 
is c orrelate d to the amount of Cepheid mass discrepancy. 

iKelled ([2008) argued that to account for a 20% mass dis- 
crepancy, stellar evolution models required an average amount 
of convective core overshooting, - 0.4 (A^ = 0.8 for his 
model s and ranging from A^ = 0.1 to A^ = 1.5). 'Neils on et al.l 
(l20Tll) found that a lower amount of convective core over- 
shooting is required, ~ 0.3, when pulsation-enhanced mass 
loss during the Cepheid stage of evolution is included in 
stellar evolut i on mo dels. This result still holds even though 
iNeilson et al.l (201 1) misinterpreted the overshooting parame- 
terization from Keller (2008), because of the large range for 
the measu red mass discrepancy. F or the case of OGLE-LMC- 
CEP0227, ICassisi & SalarisI (1201 ih fit the observed parameters 




7000 6500 6000 5500 5000 4500 4000 



Teff (K) 



Fig. 1. Stellar evolution tracks computed using the iHeger etakl 
(l2000l) stellar evolution program. The red solid line represents 
ffc = 0.2, green dashed line ac - 0.3, blue small-dashed line 
Oc = 0.4 and the magenta dotted line is Oc - 0.45. The black tri- 
angles with error bars represent the observed radius and effective 
temperature of the binary components. 

with stellar evolution models assuming an overshooting param- 
eter Oc = 0.2, which is consistent with a mass discrepancy of 
10%, based on the given definition in the Introduction. 

We compute stellar evolution models with varying values of 
metallicity and a^- to fit the physical parameters describing the 
two components of the Cepheid binary. It must be noted that 
the metallic ity and helium abundance in stellar models from the 
iHeger etal] (|2000) cod e are correlated. From th ese fits, we de- 
termine how robust the ICassisi & SalarisI (1201 ll) results are and 
if the binary Cepheid might have a mass discrepancy consistent 
with the LMC Cepheids studied b v . Keller & Wood (.2006,) . 

3. Cepheid on the first crossing 

Stars evolve across the Cepheid instability strip multiple 
times; the first crossing occurs when a star evolves along the 
Hertzsprung Gap to become a red giant star. This first cross- 
ing is usually not considered for most studies of Cepheids (i.e. 
Bono et al. 2000) because the first-crossing timescale is about 
one-tenth the blue-loop evolution timescale. Thus, it is a reason- 
able assumption that most Cepheids are evolving along the blue 
loop; however, there is only one known ecli psing binary Cepheid 
pulsating in the fundamental mode (though lLepischak & Welch! 
(20Q'i) reported the discovery of a eclipsing binary Cepheid pul- 
sating in the first overtone). Furthermore, the two stellar compo- 
nents of OGLE-LMC-CEP0227 appear to be at a similar stage 
of stellar evolution according to the masses and effective tem- 
peratures. Also, because the mass ratio is ^ = 1.00 ± 0.01, one 
cannot rule out the possibility that both stars might be evolving 
along the Hertzsprung gap. 

We compute stellar evolution tracks for stars with initial 
mass M - 4.21 Mq with varying amounts of convective core 
overshooting assumed, = 0.2,0.3,0.4, and 0.45. An initial 
mass of 4.21 Mq is chosen so that during post main sequence 
evolution the model will have a mass of 4.17 - 4.20 Mq due 
to mass loss. Furthermore, this is the largest value of the stel- 
lar mass, consistent with observations, allowing us to predict the 
minimum amount of convective core overshooting required to 
agree with the measured stellar parameters. In Fig [Tj we plot 



Neilson et al.: Cepheid binary and mass discrepancy 



3 



evolutionary tracks on the Hertzsprung-Russell diagram, along 
with the radius and effective temperature for each star in the bi- 
nary system, where the hotter star is the Cepheid. The stellar 
models evolve along the Hertzsprung gap with mass w 4.19 Mq 
due to mass loss during main sequence evolution. The stellar 
evolution tracks suggest that the radii and effective temperatures 
of the stars in the binary system are consistent with stars evolv- 
ing across the Hertzsprung gap and a first-crossing Cepheid for 
Qc = 0.4-0.5. 

While it is possible for the two stars to be evolving along 
the Hertzsprung gap, it is highly unlikely. For instance, we com- 
puted stellar evolution models for initial masses M\ - 4.22 Mq 
and Mo - 4.21 with ac - 0.4 Both evolutionary tracks 
predict effective temperatures and radii consistent with observed 
values, but the timescales do not agree. The more massive evolu- 
tion model evolves to a temperature of Teff - 5080 + 270 K from 
the zero age main sequence in about 162 Myr, while the less 
massive model has an effective temperature Tgff = 5900 + 250 K 
in 161 Myr. However, this timescale difference depend on the 
difference between the two stellar masses, hence the mass ra- 
tio q. Two models with initial masses 4.21 and 4.22 Mq has a 
value of ^ = 0.998, and reducing this mass difference will in- 
crease the mass ratio q ^ \ and reduce the timescale differ- 
ences, all consistent with the measured properties of the binary 
system. We can show this in another way, the primary star in 
the binary system is more evolved than the secondary, meaning 
that the primar y spent less time evolving on the main sequence, 
i.e. T (X M^^ dKi ppenhahn & Weigertlll994h . Furthermore, the 
difference in timescale for the secondary to evolve to the same 
location as the primary is At/t - 2.5AM/M. The difference in 
time for the Cepheid to evolve along the Hertzsprung gap to the 
same location as the primary is of the order 10^ yr and the main 
sequence lifetime is t = 1.60 Myr, hence AM/M - 0.0025 or 
q = 0.9975. Therefore, for the two stars to be evolving along the 
Hertzsprung gap with fundamental parameters consistent with 
the observations for the same age then the mass ratio q > 0.997. 
Only more precise mass determinations of the two stars such that 
q is sufficiently less than unity will rule out this scenario. 

4. Cepheid evolution with varying metallicity 

Both iPietrzvriski et al.l (1201 Oh and ICassisi & SalarisI (1201 Ih as- 
sume that the system OGLE-LMC-CEP0227 has a metallic- 
ity of Z = 0.008. Keller & Woo d (2006.) found, however, 
LMC Cephei ds have metalhcit i es ranging fr om Z - 0.005 to 
Z = 0.015. iRomaniello et al.1 (l2005ll2008l) also measured a 
similar, significant metallicity spread in LMC Cepheids. These 
works suggest that it is possible that the binary Cepheid may 
not have the standard Large Magellanic Cloud metallicity. 
ICassisi & SalarisI (|201 ll) did test stellar evolution models with 
varying metalHcity fromZ - 0.006 toZ = 0.01 with dY/dZ-lA 
and found that the change of metallicity did not affect the amount 
of convective core overshooting necessary to fit the observed 
radii. 

A secondary consideration from Fig. [T] is that the non- 
Cepheid component has a radius and effective temperature that 
is consistent with a stellar evolution model with initial mass 
4.21 Mq and a^- = 0.3. The evolution track, however, does not 
predict a large enough blue loop to compare to the radius and ef- 
fective temperature of the Cepheid. It is unclear what determines 
the width of a blue loop and is very sensitive to the initial param- 
eters of the stellar model, such as convective core oversh ooting, 
mass-loss prescription (Neilson & Lester 2008; Neils on et al.l 
1201 ll) . and helium abundance (iFiorentino et al.ii2002) . Even the 
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Fig. 2. Stellar evolution tracks computed using the iHeger etal] 
(I2OOOI) stellar evolution program with mass M - 4.21 Mq, con- 
vective core overshooting with ac - 0.2 and varying values of 
the he lium abundance, Y, and metallicity, Z. The Heger et al] 
dlOOOl) code assumes the values of Y and Z to be correlated, 
such that dY/dZ - 2. The red solid line represents Y,Z = 
0.250, 0.005, green dashed line 0.256, 0.008, blue small-dashed 
line 0.264, 0.012 and the magenta dotted line is 0.270, 0.015. The 
black triangles with error bars represent the observed radius and 
effective temperature of the binary components. 



heliu m nuclear generatio n rates affect the blue loop evolution of 
a star (iMorel et al.ll2010l) . In this respect, if we extrapolate the 
width of the blue loop for the ac - 0.3 evolutionary track to 
smaller radiii, and hotter effective temperature, then we find that 
the Cepheid is potentially consistent with ac - 0.3 as well as 
ffc = 0.2. Therefore, the binary components may be consistent 
with blue loop evolution with ac > 0.2 for Z = 0.008. This 
result is limited by the large uncertainties of the measured ef- 
fective temperatures; more precise measurements of the Cepheid 
effective temperature would constrain the minimum extent of the 
Cepheid blue loop and better constrain the input physics for stel- 
lar evolution models. 

In Fig. |2] we plot stellar evolution tracks with ac = 0.2, 
M - 4.21 Mq with varying metallicity, Z, and helium abun- 
dance, Y. It is clear that stellar evolution models predict less lu- 
min ous blue loop s for decreasing helium and metal abundances. 
Vall e et al.l (l2009l) found that the blue loop luminosity increased 
with increased helium abundance for a constant metallicity. They 
also found that varying the metallicity for a constant helium 
abundance changed the extent of the blue loop as well as led 
to a brighter blue loop for decreasing metallicity. Our models 
are consistent with this predicted blue loop behavior 

The M - 4.21 Mq, ac = 0.2 model has a temperature and 
radius consistent with the measured values for the Cepheid at an 
age of T = 150 Myr, while the same model is consistent with 
the red giant star at an age of 154 Myr. Therefore, At/t = 0.026 
hence q x 0.9889, marginally consistent with the observed mass 
ratio q - 1 .00 + 0.01 . Cases where the Cepheid mass is assumed 
to be smaller will have a greater age for the given Cepheid pa- 
rameters, likewise models with more convective core overshoot- 
ing will also have longer main sequence lifetimes and hence 
larger ages. Because the ages are greater then the predicted mass 
ratios for models that agree with the observed parameters will 
approach unity and agree with the observed q. 
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Fig. 3. Stellar evolution tracks computed using the iHeger etal] 
stellar evolution program with mass M - 4.14 Mq, con- 
vective core overshooting with Uc = 0.3 and varying values of 
the helium abundance, Y, and metallicity, Z. The line types and 
colors have the same values of Y and Z as those in Fig.|2] 

We repeat the analysis with stellar evolution models with 
mass M = 4.14 Mq and ac - 0.3 in Fig.|3] Again, we find that 
the stellar evolution tracks have blue loops that are too small for 
fit the observed parameters of the Cepheid. It is expected that 
a smaller mass stellar evolution model will have a smaller blue 
loop luminosity than a more massive model for the same amount 
of convective core overshooting. On the other hand, increasing 
the amount of overshooting in the lower mass model reduces 
this luminosity difference between the two models. Therefore, 
it appears that stellar evolution models with Oc = 0.3 and mass 
M - 4.14 Mq can fit the observations, assuming standard LMC 
metallicity and extrapolating the width of the blue loop like be- 
fore. The observed mass uncertainties of the two components 
of the binary system, while small, still allows for a variation of 
convective core overshooting. 

Furthermore, if the stars are relatively metal-rich, then the 
observed parameters can be fit with models assuming even larger 
amounts of convective core overshooting, i.e. smaller masses 
and smaller metallicities require greater amounts of overshoot- 
ing to reproduce observed luminosities. Clearly, a metallicity 
variation consistent with the observed spread of LMC Cepheids 
suggests that the amount of core overshooting necessary for stel- 
lar evolution models to fit the masses of the binary system also 
varies. Also, noting that a value of =0.1 means a mass dis- 
crepancy of 5% then this suggests that the binary system is con- 
sistent with a Cepheid mass discrepancy of 10-20%. 

5. Discussion 

We have presented two alternative scenarios to fit the observed 
parameters of the Cepheid binary OGLE-LMC-CEP0227. The 
first scenario assumes the two stars are evolving across 
Hertzsprung gap, meaning that the Cepheid is on the first cross- 
ing of the instability strip. The second scenario considers how 
fitting the observed radius and effective temperature depends on 
the assumed metallicity. 

If the stars are evolving across the Hertzsprung gap then the 
observed masses are fit only if - 0.4-0.5, hence the mass 
discrepancy would be 20-26%; assuming convective core over- 
shooting is the only solution to the mass discrepancy. If this is 



the correct scenario describing the evolution of the two stars then 
one constrains the underlying physics explaining the Cepheid 
mass discrepancy. For instance, observations of low-mass bi- 
nary stars sugg est that the convective core overshooting parame - 
ter is Oc = 0.2 (ISandberg Lacv et al.ll20Tol; IClausen et al.ll2010l) . 
leaving about 10 - 60% mass discrepancy unexplained. Because 
the Cepheid is on the first crossing pulsation-driven mass loss 
(Neilson et al. 2011) cannot explain the remaining mass discrep- 
ancy. The star has not spent enough time as a Cepheid to lose 
enough m ass. Thus rotational mixing remains a viable option 
(iBono et al. 2006; Mevnet 2009). However, we note that this sce- 
narios requires a mass ratio q > 0.997. 

The second scenario assumes that the stars are evolving 
along the blue loop with non-standard LMC metallicity and he- 
lium abundance. We find that varying the composition plays an 
important role in predicting the mass of the stars from evolution 
models, and a change in Y and Z, consistent with the spread 
of metallicities determined for LMC Cepheids, suggests that 
the amount of convective core overshooting requir e d to m atch 
observations is - 0.2-0.35. Noting that iKelled ( |2008|) and 
Neilson et al. (2011) argue that a value of Uc = 0.1 (A^ = 0.2) 
accounts for a mass discrepancy of 5% and if we extrapolate the 
width of the blue loops, then we have a measure of the mass dis- 
crepancy. In this case, we find a Cepheid mass discrepancy of 
10-18%. 

The evolutionary state of these stars may be understood in 
greater detail using measurements of the rate of period change, 
or the ni trogen abundance and by the metallicity in general. 
iTurner et al. (2006) showed that first-crossing Cepheids have pe- 
riods that are increasing with time and that the rates of period 
change tend to be significantly larger than for Cepheids on the 
third crossing of the instability strip. Measuring the nitrogen 
abundance will also provide information on whether the stars 
have evolved throu gh the first dredge-up on the red giant branch 
(Hu nter et all ^2008). For instance, if a M = 4.21 Mq LMC 
Cepheid is evolving on the first crossing then [N/Fe] 0.3 while 
evolving on the blue loop will have [N/Fe] x 0.8. It must be 
noted that rotational mixing during main sequence evolution also 
enhances the nitrogen abundance, implying that first-crossing 
Cepheids may also have nitrogen enhancements. However, ob- 
servations of both the rate of period change and nitrogen abun- 
dance would constrain the main sequence rotatio n and mix- 
ing of these stars. Similarlv. iKovtvukh et al.l (|2005') argued that 
Cepheids evolving along the blue loop will be lithium depleted 
and suggested this as another method to constrain the evolution 
status of Cepheids. The second scenario can be constrained by 
spectroscopic observations of the binary stars to determine the 
metallicity. The helium abundance is much more difficul t to de- 
termin e, helium lines disappear in cool star spectra. Dupre eTt al.l 
measured helium abundance variations in globular clus- 
ter stars based on the equivalent width of the He I /110830 
line, w hich has been observed in Cepheids dSasselov & Lesterl 
ll994ab). 

It should be noted thatl Pietrzvnski et"al] (1201 ll) reported the 
analysis of a second LMC eclipsing binary with a Cepheid 
component, OGLE-LMC-CEP-1812. This new Cepheid binary 
seems to provide another stringent test for binary star evolution, 
as its mass ratio and evolutionary phases are inconsistent with its 
age. The secondary star appears to be about 100 Myr older than 
the Cepheid primary. The Cepheid mass here is also consistent 
with pulsation masses, but no evolutionary models were tested 
in that work. 

The results presented in this work indicates that there is not 
enough information known about the binary Cepheid OGLE- 
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LMC-CEP0227 to precisely constrain stellar evolution models 
and hence resolve the Cepheid mass discrepancy. The predicted 
mass discrepancy does not rule out the Cassisi & Salaris (20ll|) 
results, but suggests that other scenarios are possible, such as 
more convective core ov ershooting and pulsation-driven mass 
loss dNeilson et al.ll20in) . Only more detailed observations will 
constrain the underlying physics that solve the di screpancy and 
test th e results of .Cassisi & Salaris (201 lb and iNeilson et alj 
(l20Tlh . 
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